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LIST  OF  SYMBOLS  * 


a  spreading  rate  parameter  =  6/x 

c  specific  heat  at  constant  pressu.e 
P 

ca  local  mass  fraction  of  the  region  1  gas 

Cp  local  mass  fraction  of  the  region  2  gas 

e  total  energy  transfer  rate  per  unit  area 

h  local  enthalpy 

6  mixing  length  (see  Eq.  3.22) 

M  local  Mach  number 

m  local  molecular  weight 

n  local  particle  density 

p  local  pressure 

q  local  heat  transfer  rate  per  unit  area 
R  local  gas  constant 

ft  universal  gas  constant 

St  Stanton  number  ~  -  h2) 

T  local  temperature 

time 

u  local  velocity  in  the  r  direction 

v  local  velocity  in  the  y  direction 

x  coordinate  parallel  to  the  initial  region  1  flow 
y  coordinate  normal  to  x 

7  ratio  of  specific  heats 

6  total  thickness  of  the  mixing  region  (at  a  given  x) 
q  radial  coordinate  =  y/5 

k  shear  stress  proportionality  constant  (see  Eq.  3.23) 
p  local  density 

t  turbulent  shearing  stress  (.-.ee  Eq,  3.22) 


*  Additional  symbols  introduced  in  the  appendix  are 
defined  in  the  appendix. 
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1 .  INTRODUCTION 

An  analytical  study  of  the  turbulent  mixing  of  two 
dissimilar  gases  in  a  free  shea^  layer  was  performed  by 
Aeronautical  Research  Associates  of  Princeton,  Inc.  for 
the  Hypervelocity  Kill  Mechanisms  Program.  It  was  felt 
desirable  to  develop  a  simplified  method  by  which  closed 
form  expressions  could  be  obtained  for  the  characteristics 
of  such  a  flow.  These  included  the  boundaries  of  the 
mixing  region,  the  location  of  the  dividing  streamline, 
the  .local  flow  velocities,  the  local  thermodynamic 
conditions,  and  the  mean  rate  of  energy  transport  across 
the  shear  layer. 

It  will  be  shown  that  although  the  effects  of 
compressibility  and  the  presence  of  dissimilar  gases  on 
either  side  of  the  mixing  region  raise  the  level  of 
complexity  of  the  problem,  the  approach  remains  the  same 
as  that  utilized  to  solve  the  single  gas  incompressible 
problem.  Therefore,  in  order  to  focus  attention  on  this 
approach,  the  incompressible  one -gas  flow  is  treated 
fir3t  (Section  3)>  then  the  modifications  required  to 
account  for  compressibility  are  shown  and  the  results, 
which  then  become  functions  of  Mach  number  and  enthalpy 
are  given  (Section  4),  and  finally,  the  treatment  and 
results  for  the  mixing  of  two  dissimilar  gases  are 
presented  (Section  5). 
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2.  THE  FREE  SHEAR  LAYER 

Before  proceeding  with  the  solution  of  the  incom¬ 
pressible  single  gas,  the  compressible  single  gas,  and 
finally  the  compressible  two-gas  problems,  a  more  detailed 
description  of  the  flow  pattern  under  study  and  the 
corresponding  model  utilized  throughout  the  analysis 
should  be  given.  Shown  belov;  is  a  sketch  of  a  cross¬ 
sec  ticn  oi  a  free  shear  layer  showing  the  cnco^  ng  free 
stream,  the  region  of  turbulent  mixing,  and  the  entrained 
flow.  Typical  velocity  profiles  of  the  free  stream  and 
mixing  layer  are  also  shown. 
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Sketch  of  Actual  Shear  Layer 

Several  modifications  to  the  above  picture  can  be 
introduced  which,  while  simplifying  the  analysis,  do  not, 
in  general,  materially  effect  the  solution.  For  one,  the 
boundary  layer  buildup  in  the  free  stream  prior  to  its 
arrival  at  the  leading  edge  of  the  hole  is  assumed  to 
have  little  effect  upon  the  subsequent  mixing  and  thus 
is  neglected.  For  another,  the  exact  self-similar  hori¬ 
zontal  velocity  profile  in  the  mixing  region  can  be 
approximated  by  a  properly  chosen  straight  line.  It  is 
recognized  that  analytical  expressions  for  the  actual 
profile  do  exist  (for  example,  reference  2);  however,  a 


external 
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mixing  region 


/ 


entrained  flow 
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linear  orcflle,  particularly  for  the  two-gas  compressible 
flow  problem,  allows  the  development  of  the  solution  in 
terms  of  closed  form  expressions  which  are  most  instructive 
and,  for  our  purposes ,  appear  adequately  accurate.  Foe 
method  by  which  the  proper  eauivalent  profile  was  obtained 
is  discussed  in  the  next  section. 

Several  assumptions  common  in  boundary  layer  theory 
were  also  made.  Specifically  they  were  that  the  gradients 
in  the  vertical  direction  are  much  greater  than  those  in 
the  horizontal  direction,  the  static  pressure  is  everywhere 
a  constant,  and  the  mean  flow  is  s'  ady. 

The  use  of  these  assumptions  leads  to  the  idealized 
model  for  a  shear  layer  shown  in  the  next  sketch.  A 
number  of  important  quantities  which  will  be  discussed  in 
the  development  to  be  presented  in  the  next  several  sections 
have  been  indicated. 

Tr 


Sketch  of  Model  for  Shear  Layer  Studies 
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3.  INCOMPRESSIBLE  SlNGLE-GAS  PLOW 

In  this  initial  analysis,  several  flow  parameters 
of  major  interest  will  be  derived  for  the  incompressible 
single-gas  shear  layer.  For  this  case  the  method  will 
be  most  transparent  and  experimental  data  are  available 
for  comparison.  Having  established  the  pattern  of  the 
development  with  this  simple  case,  it  will  not  be 
necessary  to  retrace  the  steps  in  detail  for  each  snose- 
quent  case.  In  the  ensuing  sections,  the  emphasis  will 
be  upon  the  modifications  required  for  treatment  of  the 
more  general  flows  and  the  subsequent  effects  upon  the 
various  flow  parameters  cf  interest. 


3.1.  Horizontal  velocity. 


If  one  assumes,  as  is  customary  in  turbulent  shear 
lay^r  theory,  that  the  mean  profiles  are  self-similar 
and  that  mixing  length  is  proportional  to  shear  layer 
breadth,  it  can  be  shown  that  the  total  spread  of  che 
layer  varies  linearly  with  downstream  distance  and  that 
the  velocity  and  hence  all  other  flow  parameters  remain 
constant  along  any  ray  emanating  from  the  origin.  Thus 
the  problem  can  be  expressed  in  terms  of  only  one  inde¬ 
pendent  variable  ^  where  r(  =  y/x  . 

Expressing  the  statements  that  the  horizontal 
velocity  is  linear  in  y  and  that  the  shear  layer  thick¬ 
ness  is  linear  in  x  as 


and 


we  can  write  in  terms  of  the  s5.ng.le  coordinate  q  =  y/6 
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from  which  it  can  be  shown  that 


+  3 

3 


rs1/a  = 


=  1 


u. 


u. 


and  t)0/a  =  - 


u. 


u-, 

JL 


(3.3) 


where  ,  t)p  ,  ,  and  f>  are  shown  in  the  final 

sketch  of  the  last  section. 

The  reference  velocity  Uq  in  turn  can  be  obtained 
by  utilizing  the  fact  chat  the  total  horizontal  momentum 
flux  in  the  mixing  region  may  be  written 

'*1 

f  pu2dy  »  p1u12y1  (3-3^) 

>2 


where  y1  and  are  the  limits  of  the  miring  region 

at  x  .  Transposing  to  the  coordinate  t)  and  nondimen- 
sionalizing,  this  may  be  written 


(3.5b) 


Introducing  the  relationships  (3.3)  and  (3.*0  for  u/u^ 
and  the  limits  r^/a  and  rj^/a  and  performing  the 
indicated  integration,  we  find 


(3.6) 


Finally,  returning  to  the  above  equations  for  the  horizontal 
velocity  and  mixing  region  limits,  we  can  now  write 

JL  ^  £  3 

U1  "  3  a 


(3.7) 
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r'l  1 


(3.8) 


n 


2 


a 


(3.9) 


3.2.  Vertical  velocity. 

The  vertical  velocity  anywhere  in  the  flow  can  be 
obtained  through  the  use  of  the  continuity  equation.  For 
incompressible  .flow,  this  equation  is 


+  $1  =  0 
dx  dy 

(3.10) 

Using  (3.7)  and  replacing  q  by  y/x  ,  we  car. 

write 

du  v_ 

dx  a  v.2 

rv 

(3.11) 

or  by  the  use  of  (3.10) 

*s. 

3y  a  ? 

(3.12) 

Integrating  this  last  expression  yields 

„  2 

JL  -  IL_  +  c 

n  Sa 

(3.13a) 

uf 

where  C  =  — ^  from  the  boundary  condition 
t,  =  .  Thus 

v  =  0  at 

t  ■  *  M  •  Ox)2] 

(3J.3b) 

or  by  substituting  (3.8) 
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v  ll 

Ux  =  nw  "  5  J  (3.13c) 


A  particular  vertical  velocity  that  may  he  of  interest 
is  the  one  that  exists  at  the  lower  edge  of  the  mixing 
region,  l.e.  the  velocity  at  which  the  gas  initially  at 
rest  is  being  entrained.  By  substituting  the  expression 
rtg/a  from  (3.9)  for  t\/&  in  the  above  equation,  we  obtain 


u 


1 


(3.14) 


3*3.  Dividing  stream! 1  no 

Of  particular  interest  in  several  future  computations 
will  be  conditions  on  the  dividing  streamline.  The 
6  ^iding  streamline;  of  the  mixing  region  is  defined  as  that 
streamline  above  which  the  mass  flow  is  equal  to  the  mass 

flow  from  the  external  stream  that  has  been  captured  by  the 
mixing  region.  Thus 


/pudy''Wi  (3.15a) 

y* 


where  the  subscript  »  will  be  used  to  denote  conditions 
on  the  dividing  streamline.  In  terms  of  the  variable  n 


%/a 


a 


(3.15b) 
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Again  utilising  (3.7)  and  (3.8)  for  u/ 

’T“tie  eXPTe06l0n  for  I*/3  ls  obtained  from  which  one 
deduces  that  ' 


a  “  J  V  Ja  ’  *o893  (3.16) 

Uote  tiiat  tue  plus  sign  of  the  square  root  is  to  be  taken 
since  %/a  must  lie  inside  r^a  =  -  2/3  . ) 

3*^*  ^££lactlon  of  the  free  stream, 

In  the  de^vations  for  the  vertical  velocity  and  the 
Vl  S  streamline,  it  was  assumed  that  the  free  stream 
velocity  contained  no  vertical  component  (v  =  o).  If, 
instead,  it  is  assumed  that  due  perhaps  to  the  presence  of 
a  shock  emanating  from  the  upstream  edge  of  the  hole  there 
s  a  vertical  velocity  component  to  the  free  stream,  then 
e  general  vertical  velocity  and  dividing  streamline 
expressions  are  modified  as  follows  : 


.9 


we  have 


J  Pudy  -  PiVi  -  pivrx 

y* 


(3.19) 


Nondimensio^alizing  and  integrating  as  before  yields 


n*  _  2 

a  "  '  3  + 


1  2  V1 

3  a  u. 


(3.20) 


An  interesting  quantity  which  can  now  be  computed  is 
the  magnitude  of  the  vertical  velocity  component  of  the 
free  stream  necessary  to  drive  the  dividing  streamline  to 
a  horizontal  position,  i.e.  =  0  .  It  is  quickly  seen 
from  the  above  expression  that 


(3.21) 


If  one  assumes  that  the  shear  layer  is  formed  as  a 
result  of  a  uniform  flow  across  an  orifice  in  an  otherwise 
closed  thin -walled  chamber,  then  as  long  as  the  dividing 
streamline  13  deflected  downward  the  mass  in  the  chamber 
will  increase  (see  sketch  below). 

free 

stream  ,.free  shear  layer  development 
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In  this  case,  to  the  order  of  accuracy  of  the  present 
shear  layer  model;  the  resulting  increased  pressure  in  the 
chamber  should  deflect  the  free  stream  upward  until  the 
pressure  is  just  such  that  a  shock  is  formed  which  deflects 
the  exterior  flow  by  an  amount  as  given  by  (3. 21). 

At  this  point  no  further  increase  in  mass  and  thus  pressure 
in  the  cavity  is  obtained  since  t)  =  0  and  thus  the  flow 
is  stabilized.  For  finite  thickness  walls,  the  deflection 
of  r>  necessary  tc  arrest  the  *  ill-up  and  the  associated 

vr 

can  be  similarly  computed. 


Up  to  this  point,  the  boundaries  of  the  mixing  region 

and  the  dividing  streamline  have  been  obtained  as  a 

percentage  of  the  nondimcnsional  spreading  rate  a  =  6/x  . 

To  determine  this  quantity  and  thus  the  magnitudes  of  T)  , 

tv  ,  and  tv  ,  it  is  necessary  to  introduce  the  turbulent 
2  * 

stress  present  in  the  shear  layer.  The  formulation  used 
for  this  study  is  one  which  assumes  that  the  turbulent 
stress  is  given  by  the  Prandtl  mixing  length  formula, 
namely 


(3.22) 


where  £  is  the  mixing  length.  Assuming  that  l  is 
constant  across  the  layer  at  any  station  and  is  proportional 


to  the  breadth  of  the  shear  layer  t  ,  we  can  write  by 

/Su  \ 

making  use  of  our  linear  velocity  profile  =  -g-  J 


(3.23) 


where  k  is  referred  t-j  as  the  shear  stress  proportionality 
constant.  For  an  actual  velocity  profile  (see  reference  2), 


the  line  of  maximum  shear  is  found  to  coincide  with  the 
dividing  streamline.  For  this  reason  the  characteristic 
shear  stress  just  derived  will  be  assumed  fro  act  on  the 
dividing  streamline. 

To  associate  a  spreading  rate  with  this  shear,  we 
first  introduce  the  momentum  equation 


Spuv  _  St 
By  By 


(3.24) 


and  integrate  this  expression  from  the  l>wer  boundary  of 
the  mixing  region  y0  to  the  dividing  streamline  y*  . 
There  results 


y* 

d  r  2  2^2  2^* 

—  j  pu'dy  ♦  p2u2  —  -  p»u,  —  ♦  P2u2v; 


V 

* 


(3.25) 


Applying  the  conditions  =  Tg  =  0  and  noting  that 
dy  /dx  =  v  /u  ,  this  expression  reduces  to 

■X  if  ir 


_d 
dx  o' 


2 

f  *(f)  p 


V 


a 


\V 


- 2 

ap1u1 


K_ 

a 


(3.26) 


Finally,  carrying  out  the  indicated  integrations  and 
differentiations  results  in 


a 

K 


8  T1 

Tt3  (3.27) 


where  rj^/a  is  in  general  given  by  (3.20).  For  the  case 
in  which  there  is  no  free  stream  deflection  (v^  =  0)  , 
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a 

m 


we  have,  from  (3.16),  r^/a  =  -  .0893  and  the  spreading 
rate  becomes 

a  =  13.6k  (3.28) 

It  appears  that  we  have  merely  substituted  one  unknown 
quantity  for  another.  However,  shear  stress  proportionality 
factors  have  been  experimentally  obtained  over  a  wide  range 
of  flow  conditions  both  at  ARAP  and  elsewhere.  For  Instance, 
in  reference  1  a  shear  stress  parameter  K  analogous  to  ;c 
was  obtained  from  the  measurement  of  free  axisymmetric 
turbulent  jets.  It  was  found  there  that  K  depends,  to 
first  order,  only  upon  a  suitably  chosen  local  Mach  number. 
The  Mach  number  used  in  reference  1  was  that  found  where 
the  local  jet  velocity  had  dropped  to  one-half  its  value  on 
the  centerline  for  the  same  axial  location,  this  location 
being  quite  close  to  the  point  where  the  stress  is  a 
maximum.  A  plot  of  K  as  a  function  of  local  Mach  number 
thus  defined  is  given  in  Figure  1.* 

To  test  the  suitability  ox’  utilizing  the  K  from 
reference  1  as  a  basis  for  the  k  parameter  developed  here, 
use  vTas  made  of  calculations  of  the  spreading  rate  of  an 
incompressible  shear  layer  found  in  reference  2  and  of 
(3.28)  for  the  incompressible  relationship  between  spreading 
rate  and  k  .  Reproduced  in  Figure  2  is  the  shear  layer 
profile  given  in  the  above  reference.  The  proper  equivalent 
linear  profile  also  shown  in  Figure  2  was  constructed  such 
that  the  momentum  lest  between  yQ  and  y1  ,  i.e.  the 
momentum  gained  from  y?  to  y^  ,  in  the  case  of  the  linear 
profile  is  equal  to  the  momentum  which  is  lost  by  the  exact 


*  The  values  at  K  shown  in  Figure  1  are  a  factor  of  4  less 
than  those  in  reference  1.  This  change  was  necessary  to 
compensate  for  the  difference  in  the  basic  definition  of  r 
used  here  and  in  the  referenced  report. 


iJUi jliifc UHI jJ ri ■  ►unirtd . . . . . ""  .  . .  . . . 
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profile  between  the  same  two  limits.  Thus  we  have  reauired 
that  the  momentum  transferred  be  the  same  for  each  profile, 
tf-il &  condition  leads  to  a  value  of  the  spreading  rate 
a  *  6/x  of  O.I85  for  the  equivalent  linear  profile. 
Peturning  to  (5.28)  we  find  that  k.  for  incompressible 
flow  is  then  .0119.  From  Figure  1  it  is  seen  that  the 
incompressible  value  of  K  is  .011?  which  agrees  quite 
favorably  with  the  above  result.  Thus  we  will  assume  in 
what  follows  that  the  evaluation  of  K  versus  M  shown  in 
Figure  1  is  valid  for  free  shear  layer  flow  in  which  M  is 
the  local  Mach  number  on  the  dividing  streamline  (the 

location  of  maximum  shear  for  incompressible  free  shear 
layer  flows). 

3'^.  Summary  of  the  incompressible  flow  results. 

Before  proceeding  to  more  general  flow  problems  and 
the  determination  of  the  energy  transferred  to  the  cavity 
associated  with  each  of  them,  it  might  be  useful  to 
summarize  the  incompressible  flow  results.  Table  3,1 
presents  the  general  form  of  the  several  expressions 
previously  derived.  Table  3.2  provides  the  specific  values 
of  important  quantities  for  the  case  of  v  ™  0  and 
«  =  ,0117.  1 


TABLE  3.1 

Incam^resoible  Shear  Layer  Flow  Expres 3  ions 


H 

a 


2 

x.  i  fjA 

U1  2  L\  a  / 

9  J  ul 

*--!♦/ 

7~TV7 

3  T  a  Uj 

a.  _  f  2  Z7!* 

K  1 3  \a  + 

.  s') .  a  r(Sf_  l  +  ,  s  ir 1 

3/  9  lW  9  J  +  if  j  j 

TABLE  3.2 

Specific  Values  of  incompressible  Plow  p^3mP(-or.a 

Vj  =  0  and  k  =  .0117 

uQ/u1  =  .667 

Vui  -  -577 

Vo/Uj  -  .0305 

\V"1  *  -0:i02 

a  ~  x  “  *^3  which  corresponds  to  an  angle  of  +10.6° 

_  y1 

\  "  IT  =  which  corresponds  to  an.  angle  of  +3.5° 

Y*  _ 

\  ~  x  ~  “01<54  which  corresponds  to  an  angle  of  -  0.9° 

_  y2  _ 

J'2  ”  x  «122  which  corresponds  to  an  angle  of  ~  7.]° 


COMPRESSIBLE  single-gas  flow 

Wien  the  Mach  number  in  the  free  stream  is  appreciable, 
the  density  can  no  longer  be  considered  constant  across  the 
layer.  In  the  previous  section,  the  integration  of  mass 
flux  and  momentum  flux  assumed  p  invariant.  Now  these 
Integrals  wist  be  reevaluated  in  terms  of  a  varying  density. 

4a,  Density  distribution 

In  h'hat  follows  we  will  assume,  as  was  done  in 
lelOcence  1  where  the  validity  of  the  assumption  is  discussed 
in  some  detail,  that  the  mean  local  enthalp3’  in  the  mixing 
region  id  related  to  the  mean  velocity  through  the  Crocco 
integral.  Thus  the  local  stagnation  enthalpy  may  be  written 
as  a  linear  function  of  the  local  velocity., 


h  =  Au  +  B 


(4a) 


Expanding  h°  in  terms  of  a  local  static  enthalpy,  this 
latter  quantity  can  oe  written  as 


h  =  Au  +  B  -  -• 


(4.2a) 


ana  after  inserting  the  boundary  condition  h  =  hx  when 
v  -  u^  and  h  =  hg  v?hen  u  -  0  ,  h  can  be  expressed  as 

h  =  \  *  (hl0  -  ha)^  -  (hi°  "  1^)  (j~)  (4.2b) 

Writing  now  the  local  density  as 


p  "  I&  “  l8r  ’  t^T 


(4.3) 


and  recalling  the  assumption  that  the  pressure  p  is 
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constant  across  the  mixing  layer,  we  obtain 


(4.4) 

where  h^/h^  is  the  stagnation  enthalpy  jjatio  across  the 
entire  la^er  and  h^/h^°  =  j  1  +  — ~  M^2  j  1  is  a  function 

of  the  free  stream  Mach  number  and  ratio  cf  specific  heats. 
Equation  (4.4)  provides  the  necessary  relationship  between 
the  density  and  velocity  for  use  in  the  above-mentioned 
Integrals.  Because  of  the  linear  relationship  between 
velocity  and  location  (3.3b  these  integrals  can  be 
obtained  in  closed  form. 

4.2.  Compressible  single-gas  flew  parameters. 

Not  all  of  the  flow  relationships  developed  in  the 
last  section  are  modified  by  c  .pressjbility .  The  deriva¬ 
tions  of  uhe  general  expressions  for  the  horizontal  and 
vertical  velocities 


a  „  ^  +  n 

vix  cl 


(3.3) 


r  .  V  2  v  2-, 

iUA  fji\ 

i  w  "  Uy  J 


v  _  aj 

U1  *  '2l' 

and  the  expression  for  the  mixing  .region  limits 


(3.13b) 


\  z  ^0 
a 


u 


ana 


0 


a 


u. 


(3.4) 


did  r.Oo  involve  the  density  and  thus  are  unaltered. 


However,  the  reference  velocity  uQ/al  ,  the  spreading 
rate  a/c  ,  and  the  dividing  streamline  r^/u  obtained 

"s  a  result  of  momentum  or  mass  flux  Integration  w'zt  be 
reevaluated. 


Reference  ve/oclt^.  The  velocity  uQ/u1  was  obtained  by 
equating  the  total  momentum  in  the  shear  layer  at  any 
point  to  the  captured  momentum  of  the  initial  flow  (see 
section  3.1/.  For  the  compressible  shear  layer,  the 
integral  of  equation  (3.5b)  becomes 


\/* 


I 


d 


where  the  density  ration  is  given  by  equation  (4.4). 
Integration  now  results  in  an  explicit  but  rather  involved 
expression  for  Uq/^  .  Written  functionally 


un  v 

\  ’  M1  '  yl)  <4-6) 

The  actual  expression  is  given  In  the  appendix.  It  has 
been  evaluated  for  several  enthalpy  ratios  over  a  range 
of  Mach  numbers  on  ARAF's  high  speed  digital  computer. 

The  results  can  be  seen  in  Figure  3.  Using  this  Infor¬ 
mation,  r(i/a  and  ri2/a  are  immediately  found  from  (3.4;. 

.  The  dividing  streamline  was  obtained 
Irom  mass  flow  considerations  (see  section  S.’-).  The 
integral  In  (3.15b)  when  the  density  Is  included,  yields 
an  equation  for  r|#/a  in  terms  of  the  compressible 
parameters  h2/n^  ,  mi  ,  and  7.,  and  the  vertical 


jkrj  EU'  3 1  v  r_-. : 


component  of  the  free  stream  v]/u1  .  Thus 

/Hi  h2  „  vi\ 

V  ’  ^0  ’  M1  <  Vi  »  -  0 
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(4.7) 


(Again,  the  actual  expression  is  available  in  the  appendix  j 

However  before  evaluation  of  „*  itself  (as  well  as  n, 

and  V'  the  spreading  rate  a  must  also  be  recomputed'.' 

^a^rate.  The  rate  at  which  the  mixing  region  spreads 
a  -  6/a  was  o Drained  for  the  incompressible  case  by 
computing  the  ratio  of  a  to  the  shear  stress  factor  * 

From  this  ratio  and  a  knowledge  of  *  ,  the  nondimenslonal 
sprea  .  ng  rate  6/x  can  be  obtained.  For  compressible 
flows  the  integration  Indicated  in  (3.26)  yields 


a  h 

*-hL 


M, 


*1 


a ) 


(4.8) 


(see  appendix  for  actual  expression.)  The  proper  value  of 

IT  1  O  **  _ t  t  n  . 


as  a  function  of 


" ne<3  by  first  computing  M„ 

2^1  ’  >  anu  >2  ’  anrl  then  using  the  assumed  rela- 

onship  between  k  and  ,  given  in  Figure  1.  The  local 

value  of  Mach  number  M*  is  shown  in  Figure  4.  The 
resulting  evaluation  of  a  =  6/x  is  presented  in  Figure  5. 

-MX.*?  S6en  thEt  enthalpy  r’atl°  across  the  shear  layer  has 
e  effect  on  the  total  spreading  rate, 

With  the  compressible  spreading  rate  known,  we  can 
now  return  to  (4.7)  and  (3.4)  and  obtain  the  dividing 
at*>e*ml*ne  and  mixlnS  region  limits,  respectively.  These 

and 


quantities  are  shown  as  fractions  of  Mach  number  M 
enthalpy  ratio  h^0  m  Figure  6.  1 


2.4 


4S  ihitm'ihrtfttwwwrit. 


4.3.  Energy  transfer. 
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A  characteristic  of  free  3hear  layer  flow  that  is  of 
interest  when  there  exists  a  difference  between  the  level 
of  stagnation  enthalpy  on  each  side  of  the  layer  is  the 
resulting  rate  at  which  energy  is  transferred  across  the 
layer.  To  compute  its  magnitude,  we  start  by  writing  that 
the  totaj.  energy  transferred  across  a  region  of  the  shear 
layer  per  unit  time  per  unit  area  is 


e  =  q  +  tu 

where  q  ,  the  local  heat  transfer  rate,  is 


q  =  pi2 

oy  I  cy 


and  from  our  previous  discussion  of  shear 

T  =  p'2 1  !|l 

Thus  the  total  energy  flux  can  be  written 


(4.9) 


(4.10) 


(3.22) 


e  =  p£‘ 


du 

Sy 


f  dh  3u 
t-  +  u  ™ 
Ldy  ayj 


du | oh° 
dyl  3y 


(4.11) 


Now,  utilizing  our  previous  assumptions  of  a  linear 
variation  of  stagnation  enthalpy  with  velocity  and  a 
linear  velocity  profile 


0 

top  m  du  _  (  hl  •hgyuA 

3y  du  dy  V  JK~T) 


(4.12) 


rftffiSia 


and  since  from  previously  obtained  results 
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21  aui 

dy! 


k6m., 


the  equation  for  total  energy  transfer  becomes 

*  =  KP*\  (h®  ~  h2J 


(*.13) 


Here  we  have  assumed,  as  was  the  case  for  the  shear,  that 
the  characteristic  value  of  a  given  transport  should  he 
computed  on  the  dividing  streamline. 

There  are  several  ways  in  which  this  quantity  can  be 
presented.  Commonly  the  Stanton  number  is  utilized  Here 
from  a  comparison  of  the  definition  of  Stanton  number 


e  = 


=  PjU-^t  (  h~  -  h, 


(4.14) 


with  (4.13),  it  can  be  seen  that 


P* 

St  =  JtK 
?1 


(4.15) 


The  density  ratio  on  the  dividing  streamline  may  be 
computed  by  first  calculating  u^  and  q/a"  from  (4.6) 

-  !  '"I5"  U#/U--  fr™  (3m3>’  ana  finally  p,/Pl 

*r°m  •  •  Tne  shea;  stress  proportionality  constant  k 

SlVen  ±n  the  prevloue  section.  The  Stanton  number  thus 
computed  for  the  range  of  Mach  numbers  and  enthalpy  ratios 
previously  shown  is  given  in  Figure  7, 

A  second  and  perhaps  here  a  more  useful  nondlmenslonal 
Parameter  for  the  calculation  of  energy  transfer  may  be 


kl 


obtained  a3  follows.  Expressing  the  density  p,  us 
(see  (4.3) ) 


-  1  h. 


Vl 


the  relation  for  energy  transfer  rate  (4.14)  can  be 
expressed  as 


(4.16) 


e  ~ 


*1 


u,  St 


7i  “  1 


-  ]  ul 


h2Yhl\ 

■  MV 


(<*.17) 


.Dividing  (4,17)  by  (7^/(7^  -  - } j  P, ,  one  obtain? 


tT^T  <Y‘i 


/ 

St  (1 


h. 


V'hl' 


(4.16a) 


P* 

cr  since  St  ■-  — -  k  there  results  finally 
P1 


!•(■-  3Xg) 


yprr 


(4. 18b) 


A  plot  of  e/fy^/Cy^  -  ljjp^u^  is  given  in  Figure  8. 


rr~^aegygKg¥)sfs>*--.?g¥-. 
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5.  COMPRESSIBLE  TWO -GAS  PLOW 

In  this  section  we  will  extend  the  analysis  to  include 
the  mixing  of  two  dissimilar  gases. 

5.1.  Density  distribution. 

In  what  follows,  we  again  utilize  the  development  in 
reference  1.  There  it  was  shorn  that  for  the  mixing  of 
two  species  (in  which  no  chemical  reaction  takes  place) 
the  local  mas3  fraction  of  each  species  varies  linearly 
across  the  mixlnto  region.  Letting  ca  and  c^  represent 
the  mass  fraction  of  'he  two  species  defined  as 

cc  pc/p  and  =  pp/p  (5.1) 

then  from  reference  1 


c- 


a 


u 


(3.2) 


c 


u 


(5,3) 


’where  a  represents  the  species  in  the  free  itream. 

tfith  these  relationships,  we  are  in  a  position  to 
reevaluate  the  local  density  expression.  Assuming  that  ^he 
gases  are  still  perfect 


Pcpm 


fth 


(5.4) 


where  ft  is  the  universal  gas  constant  and  m  the  local 
molecular  weight.  Thus,  for  use  in  our  mass  flux  and 


momentum  flux  integral, 


JL  _  icj>/oDl )  ( ni/rni ) 
Pi  (V^) 


I't;  n't 


from  which  we  see  that  in  addition  to  enthalpy  ratio  which 
s  already  available,  we  need  expressions  for  the  specific 
heat  and  molecular  weight  ratios. 

The  local  enthalpy  can  be  expressed  as 


C  m  =  h  —  r>  H 

P"  an 


a*a  +  cphp  -  caV,T  +  c6cd_T 


(5.6) 


if  Ideal  gases  are  assumed  as  before.  Eliminating  the 

temperature 


CP  ~  cacPjL  +  cpcp. 


or  using  (5.2)  and  (3.3) 


(5.7) 


■  ©  \ 


*  (l  -  s'.' 

\  uw 


ana  so 


Ji-A  +% 

>1  Ut  \  cn,/  C„, 


(5.8) 


(5.9) 


The  molecular  weight  can  be  expressed  as 


71  ^  £ 
n 


where  n  is  the  particle  density.  Noting  that 


pa  naml 


and  c 


3  "  p 


£  =  V!s 


(5.10) 


(5.11)  - 
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so  that 

/c  C~s 

n  -  n  +  nQ  =  (~  +  —S')  (e;  12) 

a  P  yr^  ra2y 

we  can  write  by  combining  (5.10)  and  (5.12) 


Again,  using  (5.2)  and  (5.3) 


(5.13) 


(5,14) 


(5.15) 


Returning  to  (5.5)  and  substituting  for  the  specific 
heat,  molecular  weight,  and  enthalpy,  expressions  (5.9), 

(5.15) >  and  (4.2b),  we  obtain  the  desired  expression  for 
the  density. 


(5.16) 
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5.2.  Compressible  two-gas  flow  parameters . 

in  section  4,.  it  was  demonstrated  that  once  the 
proper  density  expression  had  been  formulated  in  terms  of 
the  local  velocity,  the  expression  for  ail  of  the  flow 
parameters  of  interest  could  be  obtained  in  closed  form. 
Having  reevaluated  the  local  density  expression  to  account 
for  the  mixing  of  dissimilar  gases,  we  may  proceed  in  a 
manner  identical  to  that  outlined  in  section  4.  The 
required  integrations  are  considerably  more  involved  but 
still  tractable.  The  final  expressions,  now  functions  of 
cp0/cp,  and  mT/m2  as  well  as  h2/h®  ,  M,  ,  71  ,  and 
v^/u-j_  are  given  in  the  appendix.  These  also  have  been 
programmed  for  ARAl-s  digital  computer.  For  maximum 
molecular  weight  effects,  solutions  were  obtained  for  a 
ratio  of  entrained  gas  weight  to  free  stream  gas  weight 
of  both  one-tenth  and  ten.  In  each  case  the  entrained 
gas  molecules  were  assumed  to  possess  a  large  number  of 
degrees  of  freedom  so  thao  a  72  of  1.2  was  used 
throughout.  To  Isolate  the  effects  of  the  mixing  of 
dissimilar  gases,  all  computations  assumed  a  stagnation 
enthalpy  ratio  of  0.5. 

The  same  flow  parameters  as  were  presented  for  the 
single-gas  compressible  solut. .ns  are  shown  for  two-gas 
flow  in  Figures  9  through  14.  Note  the  very  small 
variation  of  a  =  5/x  with  molecular  weight.  The  final 
plot,  Figure  14,  indicates  that  for  the  3ame  external 
conditions,  at  high  Mach  nur,  ers,  the  presence  of  a  heavy 
entrained  gas  can  be  expected  to  decrease  the  energy 
transferred  by  approximately  65  percent  relative  to  the 
case  in  which  the  entrained  gas  is  a  factor  of  100  lighter. 


. . . mu  iiiiniiiiiwnniwimiiL.jifll 


Figure  12.  Variation  of  mixing  region  boundaries  and  dividing 

streamline  -  dissimilar  gases. 


Variation  of  Stanton  number  -  dissimilar  gases 


0.0  i  2 
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6.  CONCLUSIONS 

By  the  use  of  several  simplifying  assumptions,  closed 
form  expressions  have  been  obtained  for  the  properties  of 
a  turbulent  free  shear  layer  In  a  compressible  nonisoener- 
getic  flow  in  which  dissimilar  gases  are  present  on  the 
two  s^des  of  the  mixing  region.  Similar  expressions  are 
obtained  for  both  a  compressible  anu  an  incompressible 
single-gas  shear  layer. 

Among  the  properties  in  the  shear  layer  for  which 
expressions  have  been  developed  are 

a.  the  local  velocities 

b.  the  local  thermodynamic  quantities 

c.  the  boundaries  of  the  shear  layer 

d.  the  dividing  streamline 

e.  the  rate  at  which  energy  is  transferred  across 
the  layer. 

The  specific  equations  are  available  In  the  appendix. 
Plots  of  these  quantities  as  a  function  of  Mach  number, 
enthalpy  ratio,  and  molecular  weight  are  given  in  the 
body  ox"  this  report. 
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APPENDIX 

The  relationships  for  the  several  properties  of  the 
sLigie  gas  compressible  Bh^ar  layer  are  as  follows’ 
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For  the  dissimilar  gas  compressible  problem,  these 
expressions  are  as  follows: 
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